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1. Introduction 

The origin of the nuclear force is one of the major unsolved problems in particle and nuclear 
physics. To describe the elastic nucleon-nucleon (NN) scattering at low-energies and the deuteron 
properties, the notion of the NN potential turns out to be very useful [jl]]. The phenomenological 
NN potentials which can fit the NN data precisely are known to have the following properties: 
(i) The long range part (the relative distance r > 2 fm) is dominated by the one-pion exchange 
introduced by Yukawa (ii) The medium range part (1 fm < r < 2 fm) receives significant 
contributions from the exchange of two-pions {nn) and heavy mesons (p, ft), and a), (iii) The 
short range part (r < 1 fm) is best described by a strong repulsive core as introduced by Jastrow 
[^]. (iv) A strong attractive spin-orbit force in the isospin 1 channel exists at medium and short 
distances, (i) is related to the tensor force which is a key for the deuteron binding, (ii) is important 
for the binding of nuclei with more than 2 nucleons, (iii) is important for the stability of nuclei and 
neutron stars, and (iv) is related to the neutron pairing which leads to the neutron superfluidity 
inside neutron stars 

A repulsive core surrounded by an attractive well as seen in the phenomenological nuclear 
force is a common feature of the "effective" potentials between composite particles. The Lenard- 
Jones potential between neutral atoms or molecules is a well-known example in atomic physics. 
The potential between ^He nuclei is a typical example in nuclear physics. The origin of the re- 
pulsive cores in these examples are known to be the Pauli exclusion among electrons or among 
nucleons. The same idea, however, is not applicable to the NN potential, because the quark has not 
only spin and flavor but also color which allows six quarks occupy the same state without violating 
the Pauli principle. To account for the repulsive core of the NN force, therefore, various different 
ideas have been proposed so far [Q] : an exchange of the neutral (O meson as proposed by Nambu 
[^], exchanges of non-linear pion field, a combination of the Pauli principle with the one-gluon- 
exchange between quarks and so on. Despite all these efforts, convincing account of the nuclear 
force has not yet been obtained. 

2. NN interactions from lattice QCD 

Under the situation mentioned above, it is highly desirable to study the NN interactions from 
the first principle lattice QCD simulations. A theoretical framework suitable for such purpose was 
first proposed by Liischer [0] : For two hadrons in a finite Lr' box, an exact relation between the 
energy spectrum in the box and the elastic scattering phase shift was derived: If the range of the 
hadron interaction R is sufficiently smaller than the size of the box R < L/2, the behavior of the 
Bethe-Salpeter (BS) wave function ^f{r) in the interval R<\r\< L/2 under the periodic boundary 
conditions has sufficient information to relate the phase shift and the two-particle spectrum. The 
Liischer's method bypasses the difficulty to treat the real-time scattering process on the Euclidean 
lattice. Furthermore, it utilizes the finiteness of the lattice box effectively to extract the information 
of the on-shell scattering matrix and the phase shift. This approach has been applied to the NN 
scattering lengths in 

Recently, we have proposed a closely related but an alternative approach to the NN interac- 
tions from lattice QCD |I0|]. The starting point is the same BS wave function ^Ifir): Instead of 
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looking at the wave function outside the range of the interaction, we consider the internal region 
|r| < R and define an energy-independent non-local potential U{r,r') from i/A(r) so that it obeys 
the Schrodinger type equation in a finite box. Since U (r,r') for strong interaction is localized in its 
spatial coordinates due to confinement of quarks and gluons, the potential receives finite volume 
effect only weakly in a large box. Therefore, once U is determined and is appropriately extrap- 
olated to L ^ oo, one may simply use the Schrodinger equation in the infinite space to calculate 
the scattering phase shifts and bound state spectra to compare with experimental data. Further ad- 
vantage of utilizing the potential is that it would be a smooth function of the quark masses so that 
it is relatively easy to handle. This is in sharp contrast to the the scattering length which shows 
a singular behavior around the quark mass corresponding to the formation of the NN bound state 

0- 

Since we consider the non-asymptotic region < /?) of the wave function, the resultant 
potential U and the T-matrix are off-shell. Therefore, they depend on the nucleon interpolating 
operator adopted to define the BS wave function. This is in a sense an advantage, since one can 
establish a one-to-one correspondence between the nucleon interpolating operator and the NN po- 
tential in QCD, which is not attainable in phenomenological NN potentials. It also implies that 
the NN potential on the lattice and the phenomenological NN potentials are equivalent only in the 
sense that they give the same phase shifts, so that the comparison of their spatial structures should 
be made only qualitatively. 

3. Non-local potential from the BS wave function 

Let us consider the following BS wave function for the 6-quark state with total energy W and 
the total three-momentum P = in a finite box; ^apO' = x — y,t) = {va.c\np{y,t)pa{x,t)\W,P = 
0) = Yap(f)^ '^'- The local composite operators for the proton and the neutron are denoted by 
Pa{x,t) and np{y,t) with spinor indices a and j8. The state \W,P = 0) is a QCD eigenstate with 
baryon number 2 and with the same quantum numbers as the pn system. One should keep in mind 
that \W,P = 0) is not a simple superposition of a product state |p) (g) |n), since there are complicated 
exchanges of quarks and gluons between the two composite particles. 

The spatial extent of the NN interaction in QCD is short ranged and is exponentially sup- 
pressed beyond the distance R > 2 fm. Therefore, the spatial part of the BS wave function in 
the "outer region" (r > R) satisfies the Helmholtz equation below the pion production threshold, 
(V^ + k^)Yap{^) = 0» up to an exponentially small correction. Then we can define the non-local 
potential U from y and measured on the lattice: 

{E-Ho)YapA^) = fUap-rs{7,r')YrsA7')d'r', (3.1) 
Uap-rs{r,r) = Vap-ys{f,v)S{r-f), (3.2) 

where v(= p/jj. = — /V//i) is the velocity operator. To make a formal resemblance with the 
non-relativistic case, we have introduced the "effective center of mass energy", E = /{2}Jl) = 
/mj^ and the "free Hamiltonian", Hq = — V^/m/y. By construction, the solution of Eq.(^?T]) with 
U(^p.ys{r,r') extrapolated to L ^ oo reproduces the correct BS wave function in the asymptotic 
region, and hence the phase shifts and binding energies of the two-nucleon system. 
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The simplest interpolating operators for the nucleon N = (/?,«) in terms of the quark field 
q{x) would be Np{x) = Sabc{qa{x)CYsqb{x))qcfi{x), with a, b and c being color indices and C 
being the charge conjugation matrix. Such a local operator is most convenient for relating the BS 
wave function to the four-point Green's function and the scattering observables at L — > oo. Closely 
related formulation was given long time ago by Nishijima, Zimmermann and Hagg who derived 
the generalized reduction formula for local composite fields [p^]. 

In principle, one may choose any composite operators with the same quantum numbers as the 
nucleon to define the BS wave function. Different operators give different BS wave functions and 
different NN potentials, although they lead to the same observables. This is quite analogous to the 
situation in quantum mechanics where the unitary transformation of the wave function changes the 
structure of the potential while the observables are not modified. A theoretical advantage of our 
approach based on lattice QCD is that we can unambiguously trace the one-to-one correspondence 
between the NN potential and the interpolating operator in QCD as we mentioned. 

The general form of the non-local NN potential U or equivalently the velocity dependent 



NN potential V in Eq.(|3.2|) in the two-component spinor space has been classified by Okubo and 



Marshak QISQ. The leading order (LO) and the next-leading-oder (NLO) terms of the the velocity 



expansion of V(r, v) reads [jlj] 

V = Vc{r)+VT{r)Sn + VLs{r)L-S + 0{v^), (3.3) 
= Vo(r) + Vff(r)(ai • 02) • Zi) +Vcy,{r){Oi ■ 02){ti ■ T2) 
+ [VT0{r)+VT,{r){TyT2)]Si2 

+ [VLso(r) + VL5r(r)(Ti • T2)] L • 5 + 0(v2), (3.4) 

where Vc and Vj are LO {0{^)) terms, while Vis is a NLO (0(v)) term. On the lattice, it is rela- 
tively unambiguous to extract information for the orbital angular momentum states £ = 0, 1,2,3 = 
S,P,D,F using the irreducible representations of the cubic group [Q]. Then, at most 16 independent 
(14 diagonal and 2 off-diagonal) matrix elements of the potential are obtained, so that 8 unknown 



LO and NLO terms in Eq.(3.4) can be extracted in two different ways. 



4. Central and tensor forces from lattice QCD 

To define the BS wave function on the lattice, we start from the four-point correlator. 



vac 



nfi{y,t)pa{x,t)/ (k.-f) vac) ^ Aq Wap{r-f) e-^«('-'») {t » t^), (4.1) 



where Aq is an r-independent constant. The states created by the source ^ have the conserved 
quantum numbers, {J,J^) (total angular momentum and its z-component) and P (parity). For study- 
ing the nuclear force in the = 0+ (^Sq) channel and the = 1+ (^Si and ^Di) channel, we adopt 
a wall source with the Coulomb gauge fixing at ? = Iq. The BS wave function in the orbital S -state 
is defined with the projection operator for the orbital angular momentum (P^'^') and that for the spin 

(PW) as VA(r;' So) = P^^=^^P^'=^^W^^'^) and v/-(r;3 Si) = p(^=°)p(-'=i V(^; 1+)- 

The asymptotic momentum k for the S-states is obtained by fitting the BS wave function 
i/A(r) with the Green's function G{r;k^) in a finite and periodic box satisfying (V^ + k'^)G{r;k'^) = 
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~^at('') with 5iat(r) being the periodic delta-function. The fits are performed outside the range of 
the NN interaction determined by V^VA(r)/i/A(r) JlSl]. The NN scattering lengths for the S -states 
can be deduced from the standard Liischer's formula [Q]. 

In the LO of the velocity expansion, only the central potential Vc(r) and the tensor potential 
VT{r) are relevant: The central potential acts separately on the S and D components, while the 
tensor potential provides a coupling between these two. Therefore, we consider a coupled-channel 
Schrodinger equation in the 7^ = 1+ channel [|I6|]: 

{Ho + Vc{r)+VT{r)Sy2)w{r, 1+) =£i/^(r; 1+). (4.2) 

Projections to the S-wave and D-wave components are obtained as y/„^(r;^ Si) = P^^^^^^Yapi^^'^ 1^) 
and Yapi^''^^i) = (1 ~^^^"°^)V^aj3(^' 1^)- the LO of the velocity expansion, it is sometimes 
useful to define the "effective" central potential Vf{r) [g]: Vf{r) = £ + J-^^. Note that 
y^^(r) in the ^Si channel contains the effect of Vjir) implicitely as higher order effects through 
the process such as ^Si — >^ Di Si. 

5. Numerical results in quenched QCD 

In the quenched simulations, we employ the standard plaquette gauge action on a 32^ lat- 
tice with the bare QCD coupling constant j3 = = 5.7. The corresponding lattice spacing is 
I/fl = 1.44(2) GeV (a ~ 0.137 fm) determined from the p meson mass in the chiral limit. The 
physical size of our lattice then reads L ~ 4.4 fm. We adopt the standard Wilson quark action 
with the hopping parameter k: = 0.1640,0.1665,0.1678, which correspond to nij^ ^ 731,529,380 
MeV, respectively. The periodic boundary condition is imposed on the quark fields along the spa- 
tial direction, while the Dirichlet boundary condition is imposed in the temporal direction at the 
time-slice t = 0. The wall source is placed on the time-slice at to/a = 5 with the Coulomb gauge 
fixing att = to. The lowest effective cm. energy E in the above setup ranges from —0.4 MeV to 
— 1 MeV. Note that E for scattering states can be negative in a finite box. 

5.1 Central and tensor forces in the ^Si channel 

Shown in Fig. |I| is the central potential Vc{r) and tensor potential Vt (r) together with effective 
central potential V^^{r) in the ^Si channel obtained in the LO velocity expansion. In the real 
world, V^*(r) is expected to acquire sufficient attraction from the tensor force. This is the reason 
why bound deuteron exists in the ^Si channel while the bound dineutron does not exist in the 'Sq 
channel. Now, we see from Fig. |l| that the difference between Vc{r) and V^^{r) is still small in our 
quenched simulations due to relatively large quark masses. 

The tensor potential W (r) in Fig. |l| shows that it is negative for the whole range of r with a 
minimum at short distance below 0.5 fm. If the long range part of the tensor force is dominated 
by the one-pion exchange as expected from the meson theory, Vrir) could be rather sensitive to 
the change of the quark mass. As shown in Fig. ^ it is indeed the case: Attraction of Vrir) is 
substantially enhanced as the quark mass decreases. For practical applications in nuclear physics, 
it is more useful to parametrize the lattice results by known functions. We have tried such a fit 
for Vrir) under the assumption of the one-p -exchange -i- one-pion-exchange with Gaussian form 
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Figure 1: The central potential Vc{r) and the tensor potential VV(r) obtained from the 7^ 1+ BS wave 
function at nijt = 529 MeV. 
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Figure 2: Quark mass dependence of tensor potential. The lines are the four-parameter fit using the one-p- 
exchange + one-pion-exchange with Gaussian form factors. 
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factors: V^W = b,{\ -e-'^^-^f (l + ^7 + ^ + ^3(1 -e-^^'^f {\ + ^^ + 

The results are shown by the solid lines in Fig. ' 

5.2 Velocity dependence of the potential 

So far we have considered the potential determined from the lattice data taken almost at zero 
effective cm. energy £ ~ MeV. If the local potential determined from the other energies have 
different spatial structure, it is an indication of the velocity dependent terms. Such a velocity 
dependence has been recently studied by changing the spatial boundary condition of the quark field 



from the periodic one to the anti -periodic one [17]: On a 32 x 48 lattice with the lattice spacing 



a = 0.137 fm, 2000 gauge configurations are accumulated. The minimum momentum is given by 
= (7r,7r,7r)/(32a), which leads to \pmm\ — 240 MeV and £" ~ 50 MeV. In Fig.|, the central 
NN potential for the ^Sq state with APBC ( £ ~ 50 MeV) is plotted as a function of r at f/a = 9, 
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Figure 3: The central NN potentials for the 'So state with APBC (red bars) and PBC (blue crosses) in 
quenched QCD at f/a = 9. 

together with the one with PBC {E ~ 0). Fluctuations of the data with APBC at large distances 
(r > 1.5 fm) are mainly caused by contaminations from excited states, together with statistical 
noises. The potential at r < 1.5 fm, on the other hand, is less affected by such contamination. As 
seen from Fig. ^, the NN potentials are almost identical between £ ~ and £ ~ 50 MeV. Namely, 
the non-locality of the potential with our choice of the interpolating operator is small and the LO 
potentials shown in the present paper can be used in the energy region at least up to £ 50 MeV 
without significant modifications. 

6. Numerical result in (2+l)-flavor QCD 

To compare our results with empirical data, a key role is played by a full QCD calculation on 
a large volume with a smaller quark mass. The PACS-CS collaboration is generating (2+l)-flavor 
gauge configurations by employing the Iwasaki gauge action at j3 = 1.90 on 32^ x 64 lattice and 
the 0(a)-improved Wilson quark (clover) action with a non-perturbatively improved coefficient 
csw = 1-V15 [p^]. The lattice scale is determined by niji, niK and mQ, which leads to a ~ 0.091 
fm. Thus, the spatial extension amounts to L ~ 2.90 fm. The periodic boundary condition is 
imposed along the spatial direction, while the Dirichlet boundary condition on the time-slice t = 32 
is imposed along the temporal direction. The wall source on the time-slice is located at f = with 
the Coulomb gauge fixing. Fig. ^ shows the full QCD results of the central force for rriji ~ 702 
MeV: Vcir-^So) and Vf{r,^Si) are obtained from BS wave functions on the time-slices t = S and 
t = 9, respectively, where the ground state saturations are achieved within error bars. Similar to the 
quenched results, a repulsive core surrounded by an attractive well can be seen in full QCD. 

7. Concluding remarks 

We have discussed the basic notion of the nucleon-nucleon potential and its field-theoretical 
derivation from the equal-time Bethe-Salpeter amplitude in QCD. By construction, the non-local 
potential defined through the projection of the wave function to the interaction region (the inner 



7 



Nuclear Forces from Lattice QCD 



Tetsuo Hatsuda 



3500 
3000 
2500 -, 



> 2000 



o 
> 



1500 - 
1000 - 
500 
V 



0.0 



150 
100 
50 

-50 



Vc(r;'So) 



m^=702 MeV 




0.5 



1.0 1.5 
r[fm] 



2.0 



0.0 0.5 1.0 1.5 2.0 2.5. 



2.5 



Figure 4: Central and effective central potentials in (2+l)-flavor lattice QCD simulations with niji 
MeV. 
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region) correctly reproduces the asymptotic form of the wave function in the region beyond the 
range of the nuclear force (the outer region). Thus the observables such as the phase shifts and 
the binding energies can be calculated after extrapolating the potential to the infinite volume limit. 
Non-locality of the potential can be taken into account successively by making its velocity expan- 
sion, which introduces the velocity-dependent local potentials. The leading order terms of such 
velocity expansion are the central potential and the tensor potential, and the next-to-leading order 
term is the spin-orbit potential. 

As an exploratory study, we carried out quenched lattice QCD simulations of the two-nucleon 
system in a spatial box of the size (4.4 fm)^ with the pion mass mjt = 380,529,731 MeV. The 
NN potential calculated on the lattice at low energy is found to have all the characteristic features 
expected from the empirical NN potentials obtained from the experimental NN phase shifts, namely 
the repulsive core surrounded by the attractive well for the central potential. As for the tensor 
potential obtained by the coupled channel treatment of the ^Si -state and the ^Di -state, appreciable 
attraction at long and medium distances is found. Phenomenological fit of the tensor potential 
strongly suggests the existence of the one-pion-exchange in its long range part. 

There are a number of directions to be investigated on the basis of our approach: 

1. Determination of the velocity dependence is important in deriving the NN potentials useful 



for the wide range of scattering energies. As mentioned in Sec.p.2|, studies along this line 
using the anti-periodic boundary condition in the spatial direction have been already started 



[17]. 



2. 



To derive the realistic NN potentials on the lattice, it is necessary to carry out full QCD 
simulations with dynamical quarks. As mentioned in Sec.^, studies along this line with 
the use of the (2-1-1 )-flavor QCD configurations generated by PACS-CS Collaboration are 



currently under way [16] 



3. The hyperon-nucleon (YN) and hyperon-hyperon (YY) potentials are essential for under- 
standing the properties of hyper nuclei and the hyperonic matter inside the neutron stars. 
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Recently, studies of the YN potential in quenched QCD [19| and in full QCD [2C] have been 
started. 

4. The three-nucleon force is thought to play important roles in nuclear structures and in the 
equation of state of high density matter [pT]]. Since the experimental information is scarce, 
simulations of the three nucleons on the lattice may lead to the first principle determination 
of the three-nucleon potential in the near future. 
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